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PART  ONE 

THEORY  AND  EXAMPLES 

I . INTRODUCTION 

In  a previous  report  [1]  It  was  demonstrated  that  for  a reactively- 
loaded  waveguide-fed  aperture  antenna  array  the  radiation  characteristics 
can  be  controlled  by  feeding  one  aperture  and  reactlvely  loading  the  other 
apertures.  The  reactive  loads  can  be  realized  by  placing  electrical  short 
circuits  In  the  unfed  waveguides  at  variable  distances  from  the  apertures. 
By  varying  these  short  circuit  positions,  the  antenna  beam  can  be  steered. 

The  use  of  thin  apertures  with  waveguides  of  the  same  cross  sectional 
dimensions  resulted  In  a sizable  mismatch  between  the  wave  admittance  and 
the  half-space  admittance  at  the  aperture  boundary.  This  resulted  In  small 
reactive  load  magnitudes  (referenced  to  the  wave  admittance)  and  small 
changes  in  short  circuit  distances  per  degree  change  in  antenna  beam  posi- 
tion. It  also  necessitated  the  use  of  non-standard  waveguide.  In  this 
report  a wide  aperture  antenna  array  is  analyzed.  The  above  mentioned 
problem  areas  are  reduced  by  changing  from  thin  to  wide  apertures. 

The  aperture  antenna  array  studied  consists  of  closely  spaced  aper- 
tures which  allow  large  mutual  coupling  to  exist  between  apertures.  Each 
aperture,  therefore,  has  a large  effect  on  the  resultant  antenna  beam  shape 
and  position. 

II.  STATEMENT  OF  THE  PROBLEM 

Figure  1 Illustrates  a seven  element  waveguide-fed  dlelectrlc- 
fllled  rectangular  aperture  array.  The  following  assumptions  are  made; 

1.  A perfectly  conducting  plate  covers  the  entire  z = 0 plane 
except  for  the  apertures. 

2.  The  apertures  are  one-half  wavelength  long. 

3.  The  dominant  mode  Is  assumed  to  be  the  only  mode  present  In 
each  waveguide  (all  other  modes  are  cut-off). 
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4.  The  waveguide  dielectric  filler  lowers  the  cut-off 

frequency  for  the  dominant  mode  so  that  the  assumption 
of  a cosine  tangential  electric  field  existing  across 
each  aperture  Is  reasonable. 

The  general  method  of  solution  [2]  is  to  cover  the  K apertures 
with  a perfect  conductor,  to  place  magnetic  current  sheets  +M^  and 
respectively,  on  the  left-hand  and  right-hand  sides  of  the  conductor, 
to  obtain  an  Integral  equation  for  each  (j  = 1,2,...,K)  by  equating 
the  tangential  magnetic  fields  on  both  sides  of  the  conductor,  and  to 
solve  these  Integral  equations  using  the  method  of  moments.  The  testing 
functions  are  the  same  as  the  expansion  functions  . Each  is  a 
cosine  function  In  the  direction  of  current  and  a pulse  function  in  the 
direction  perpendicular  to  the  current.  For  our  problem  the  current  is 
in  the  y direction  only. 


III.  SUMMARY  OF  BASIC  THEORY 

The  orientation  of  the  apertures  and  the  coordinate  system  is 
shown  In  Fig.  2.  The  equivalent  magnetic  current  in  the  jth  aperture 
region  is  represented  by 


= V,  M 

j 

where  the  are  constants  to  be  determined.  The  solution  given 
by  [ll  is 


(1) 


(2) 


where  V is  the  column  vector  of  the  Vj . The  magnetic  current  expansion 
function  in  the  jth  aperture  Is  taken  to  be 


Mj  = Pj  (x)  cos  ky 


(3) 


In  (2)  the  diagonal  matrix  [y'^^]  is  generated  by  two  distinct  expres- 
sions: 

K+1 

a)  For  1 = -j—  (K  is  the  number  of  apertures) , 


^11  * H^^^(cos  ky)>] 

= Y (dominant  mode  wave  admittance) 
o 


(4) 
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(Note:  < > Is  the  mathematical  symbol  for  symmetric  product 


<A,B> 


A • B dS.) 


surface 


, . „ . , K+1 

b)  For  i f -5—  , 


= [~  <-  cos  ky,  (cos  ky)>] 

Xi  WL* 

» - 1 Y cot  k,d^(0) 
o i 

where  d^(6)  is  the  distance  of  the  short  circuit  from  the 

aperture  in  the  ith  waveguide  (dependent  on  beam 
steering  angle  0) 

k is  the  free-space  wave  number 
k is  the  cut-off  wave  number 


(5) 


= /e  k^  - k^ 
r c 

is  the  relative  permittivity  of  the  dielectric  filler 
is  the  free-space  intrinsic  impedance 


Y = - /e  - (k  /k)^ 
o n r c 

In  (2)  the  matrix  [Y^^]  is  given  by 


hs 

ij 


aperture 


- Ml  • 

1 -a 


(6) 


where  (M  ) is  the  tangential  magnetic  field  in  the  ith  aperture 

generated  by  the  equivalent  magnetic  current  in  the  jth  aperture. 

In  (2)  the  column  vector  I^  consists  of  one  non-zero  element  which 
corresponds  to  the  excitation  (externally  driven  waveguide).  In 
(3)  the  expression  is  given  by 
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I: 

I 


f 

[ 

r 

i- 


i 


P^(X) 


Pl8-  2) 

0 , all  other  x. 
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The  gain  (ratio  of  the  maximum  radiation  Intensity  In  a given 
direction  to  the  maximum  radiation  produced  In  the  same  direction  from 
an  Isotropic  radiator  with  the  same  power  input)  corresponding  to  the 
direction  (0,4i)  in  the  half-space  z > 0 is  given  by  (*  means  complex 
conjugate) 

,2  ^ , 

G(e,<J>) jp”(6,4,)Vp  (7) 

8-trri  Real  (V[Y  ]*  V*) 


where 


Pj(0.<^)  = 2 


M u 

j -y 


jth 

aperture 


u 

~m 


(8) 


(Note:  P"*(0,()>)  represents  a vector  of  short  circuit  currents  when  the 

array  Is  excited  by  a plane  wave  from  the  direction  (0,<{i).  P (0,()))  can 
be  polarized  In  the  0 or  (}i  direction.)  For  the  set  of  all  P™(0,(|))  vectors 
In  the  E-plane  (0  ^ 0 £ 180°,  (j>  = 90°)  and  y-polarlzed,  we  have 


P®  (E-plane)  = p“^  (0,  f) 


/“w  J^x^cos0^sln  (y  cos  0) 


.kw 
'■  2 


Iw 

2 


CVW  A 

cos  0 


') 


(9) 


llR 

IV.  HALF-SPACE  ADMITTANCE  MATRIX  [Y  ] DETERMINATION 

In  this  section  the  application  of  a numerical  Integration  tech- 
nique to  determine  the  half-space  admittance  for  a wide  aperture  Is 
presented.  The  procedure  followed  is  to  first  subdivide  the  equivalent 
magnetic  current  sheet  over  each  aperture  region  into  p current  strips 
(see  Fig.  3).  By  using  these  current  strips  as  sources,  the  self  and 
mutual  admittances  of  these  sub-aperture  regions  are  found  by  using  pre- 
viously derived  formulas  for  the  thin  aperture  case  [1].  Finally,  these 


MAGNETIC 
CURRENT  SHEET 


MAGNETIC 
CURRENT  STRIPS 


3.  Subdivision  of  a magnetic  current  sheet  into  p equal  magnetic 
current  strips. 


yield  a single  self  admittance  term  for  the  entire  aperture 

- 

n all  of  the  aperture  of  the  array  We  o 

aperture  case  and  then  extend  the  single 

array.  -l^^Ple  aperture 

In  the  preceding  section,  the  half-space  admittance  for  a single 
aperture  was  defined  to  be  single 


where 


,hs 


M * dS 


aperture 


= iiy  V; 


HyV^  cos  ky, 


(10) 


(11) 


Let 


where 


? M 


= 1 ^ 


M 


-yV' 


W V7 

2-^-2 


(see  Fig.  3) 


^ ky 


(12) 


(13) 


- 2 + (q-1)  ^ < X < - I +(q)" 


® 1»2 


Then 


= - 2 


! 


L/2  -w/2+(q)w/p 


dy 


q-lril  i - J 

L/2  ~w/2+(q-l)w/p 

where  the  factor  2 Is  due  to  image  theory  (Th^ 

cneory.  (The  magnetic  current 

^co^lete  o»d„c„„s  .l.„,  sl„«  the  aperture  Is  covered  by  , conductor  - 
the  adMttanc.  obtained  ualng  expanaion  functions  2 M radiating 
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Into  free  space  everywhere.)  H (M  ) is  the  tangential  magnetic  field 
due  to  radiating  into  free  space.  is  approximated  by  a filament 

Uy*y^  Cp  cos  ky  of  magnetic  current  at  x = -^  + (q  - 0.5)  p By  sub- 
stituting this  approximation  for  ^ and  a similar  approximation  for 
into  (14),  we  obtain 


= - 2 (4-)  (-)^ 

wL  p 


L/2 


q=l  r=l 


-L/2 


cos  ky  • 


(15) 


qr 


where  is  the  tangential  magnetic  field  due  to  the  filament  u cos  ky 

qr 

of  magnetic  current  at  x = 4 + (r  - 0.5)  - evaluated  at  x = -^  + (q  - 0.5) 


To  solve  (15),  we  substitute  an  equivalent  expression  for  H for 

"•t 


qr 

which  the  resultant  integral  has  a known  solution.  This  expression  can  be 
derived  as  follows.  For  a source  free  region, 


VxE  = -juv]H 
V X H = jtoe  E 


Let 

and 


where 


E = - V X F 
H = V X A 


A(r) 

V- 


dr  ’ 


(16) 

(17) 

(18) 

(19) 

(20) 
(21) 


By  substituting  (18)  into  (16)  and  (19)  into  (17),  we  obtain 


»l 
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H = V X V X F 

VA  jUJV)  - ~ -A 


(22) 


F = V X V X A 

-A-  JOJVI  ~ 'A  A>- 


(23) 


Therefore,  H due  to  M is  —j  times  F due  to  J.  Using  this  relationship, 
(15)  transforms  to 


r"'. ! ! 1- 

n p q=l  r=l 


L/2 


-L/2 


dy  u cos  hy  • R 

^ AAy  ^ AAt 


(24) 


qr 


where  E is  the  tangential  electric  field  due  to  a filament  u cos  ky 


of  electric  current  at  x = - -r  + (r  - 0.5)  — evaluated  at  x = - — + 

2 p 2 

(q  - 0.5) 


For  the  self  admittance  terms  (q=r) , we  separate  the  field  and 
source  filaments  by  a distance  (~)/^  which  is  the  "equivalent  radius" 
for  a rectangular  strip.  The  quantity  in  brackets  in  (24)  represents 
the  self  and  mutual  impedances  of  thin  dipoles  in  an  array.  Therefore, 
(24)  can  be  written  as 


Y 


hs 


(-y)  (f)  (-7) 

n ' p 


\ I 

q=l  r=l 


Z 

qr 


(25) 


To  extend  this  development  to  the  multiple  aperture  case,  (10) 
can  be  rewritten  as 


aperture 


(M,)dS  . 

-^ij  -J 


(26) 


^ (i  y j)  represents  the  mutual  admittance  between  the  1th  and  jth 

of  the  1th 
c field  in 

the  ith  aperture  generated  by  the  equivalent  magnetic  curren.:  sheet  in 
the  jth  aperture. 


hs 

apertures  while  Y.  (i  = j)  represents  the  self  admittance 
^ns 

(jth)  aperture,  h"  (M  ) represents  the  tangential  magneti 

-tij  -j 
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Let 


M . = u ■—  cos  ky  , 

-i  -Iq  «-y  VV7L  ^ ' 


q=l 


(27) 


w w 


(x^  is  the  distance  from  the  origin  to  the  center  of  the  ith  aperture.) 
where 


M,  = u J-~  cos  ky, 

v-lq  ~y  V wL  ’ 


(28) 


- y + X + (q-1)  «.  < X < -^  + X + q 5^- 
2 1 ^ p—  — 2 i^p 


q — 1,2,...,  p. 


Substituting  (27)  into  (26)  we  obtain 


- - 2 
Ij 


jL  ill  I 


L -w  . . w 

T — + ^i+q  - 


2 2 
dy 


dx  • H 
^iq  v^t 


-L  -w  . , . ,\  w 

2 x^+(q-l)  - 


(iq)(jr) 


L 

2 


= (-y)  (^)  (-^)  f f ( - f dy  u cos  ky  • \ 

n p q=l  r=r  \ ^ *^(lq)(jr)  / 


2 


( 2>(l)(  2^  ^ i ^(iq)(jr) 

n p q=l  r=l  ' 'J  ' 


(29) 


where  (^  ^ J)  represents  the  mutual  impedance  between  the 

(Iq)th  and  (jr)th  dipole  current  filaments. 
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V.  REPRESENTATIVE  COMPUTATIONS 

A computer  program  using  the  preceding  formulas  has  been  written  and 
Is  Included  In  report  [1]  with  the  exception  of  YHSPA  (generates  the  half- 
space admittance  matrix  for  a wide  aperture  array  using  a numerical  Inte- 
gration technique)  which  Is  Included  In  Part  TVo  of  this  report.  For  this 
section  results  will  be  given  for  N = 7,  9,  and  15  element  aperture  arrays. 
In  all  of  the  cases  to  be  presented,  the  aperture  width  Is  0.25X  and  six 
magnetic  current  strips  are  used  per  aperture  region.  The  Initial  reac- 
tive loads  required  for  the  optimization  subroutine  are  the  reactive  loads 
which  resonate  complex  equivalent  sources.  All  of  the  gain  patterns  are 
taken  In  the  E-plane  (0  ^ 6 £ 180°,  (p  = 90°). 

Figures  4 and  5 Illustrate  the  maximum  gain  that  can  be  obtained 
for  a seven  element  aperture  array  using  a single  feed  and  reactlvely 
loading  the  other  apertures.  Tlie  triangular  symbols  A In  the  two  figures 
represent  the  maximum  gain  that  can  be  obtained  for  the  same  aperture 
antenna  by  feeding  each  aperture  with  the  optimum  complex  equivalent  volt- 
age. It  should  be  noted  that  both  excitations  yield  excessively  large  beam- 
widths  at  low  angles  (0  ^ 6 ^ 40°). 

Figure  6 represents  the  same  case  of  a seven  element  aperture  array 
considered  In  Figures  4 and  5 (for  0=0,  30°,  60°,  and  90°),  with  a IX 
loss  factor  added.  In  addition  to  the  loss  of  gain  In  the  desired  direc- 
tions, the  backlobe  In  the  6=0  pattern  Is  greatly  reduced. 

Figure  7 Illustrates  the  maximum  gain  that  can  be  obtained  for  a 
nine  element  aperture  array  using  a single  feed  and  reactlvely  loading 
the  other  apertures.  A general  Increase  In  gain  Is  observed  with  no 
apparent  change  In  the  sldelobe  structure  compared  to  the  seven  element 
aperture  array  case. 

Figure  8 Illustrates  the  maximum  gain  that  can  be  obtained  for  a 
fifteen  element  aperture  array  using  a single  feed  and  reactive  loading 
the  other  apertures,  with  a 0.1%  loss  factor  added.  Relative  sldelobe 
levels  at  6 » 60°,  75°,  and  90°  have  decreased  In  comparison  to  the  seven 
and  nine  element  aperture  arrays. 
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Fig.  A.  Radiation  gain  patterns  for  a seven  element  aperture  array.  The 
solid  lines  represent  the  pattern  for  the  reactively  loaded  array 
with  only  the  center  element  fed.  The  triangles  represent  the 
pattern  for  the  array  when  all  elements  are  fed  with  the  excitation 
for  maximimi  gain.  Aperture  width  is  0.25X.  Six  current  strips  per 
aperture  were  used  to  evaluate  the  admittance  elements. 
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. Radiation  gain  patterns  for  a seven  element  aperture  array.  The  solid 
lines  represent  the  pattern  for  the  reactively  loaded  array  with  only 
the  center  element  fed.  The  triangles  represent  the  pattern  for  the 
array  when  all  elements  are  fed  with  the  excitation  for  maximum  gain. 
Aperture  width  is  0.25X.  Six  current  strips  per  aperture  were  used  to 
evaluate  the  admittance  elements. 
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6.  Radiation  gain  patterns  for  a seven  element  reactlvely  loaded 
aperture  array  with  only  the  center  element  fed.  Aperture  width 
is  0.25X.  Six  current  strips  per  aperture  were  used  to  evaluate 
the  admittance  elements.  A loss  factor  of  0.1%  was  added  to  the 
real  parts  of  the  diagonal  admittance  elements. 
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8.  Radiation  gain  patterns  for  a fifteen  element  reactlvely 
loaded  aperture  array  with  only  the  center  element  fed. 
Aperture  width  is  0.25X.  Six  current  strips  per  aperture 
were  used  to  evaluate  the  admittance  elements.  A loss 
factor  of  0.1%  was  added  to  the  real  parts  of  the  diagonal 
admittance  elements. 


Table  1 lists  the  reactive  load  values  required  to  steer  the  beam 
from  0 to  90®  In  10®  steps  for  the  seven  element  lossless  aperture  array 
(see  Figs.  4 and  5).  (Note  that  both  Inductive  and  capacitive  reactances 
are  required.) 

Table  2 lists  the  reactive  load  values  and  equivalent  short  circuit 
distances  required  for  a nine  element  lossless  array  (see  Fig.  7).  The 
distance  magnitudes  are  referenced  to  WR-75  (2.159  cm  x 1.206  cm  — wall 
thickness  0.127  cm)  filled  with  a dielectric  = 1.55  and  operating  at  a 
frequency  of  7.88  GHz.  Admittedly,  If  larger  waveguide  and  a lower  operat- 
ing frequency  were  chosen,  the  changes  In  short  circuit  distances  required 
for  beam  steering  would  be  larger.  The  choice  of  WR-75  and  a frequency  of 
7.88  GHz  was  made  to  illustrate  that  very  small  changes  In  electrical 
short  positions  are  required  to  steer  the  beam  In  increments  of  15®  for 
waveguides  of  these  relative  dimensions  and  operating  frequencies  at  or 
above  7.88  GHz. 
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Table  1.  Reactive  load  values  (imnhos)  for  a seven  element  aperture  array 
using  six  strips  per  aperture  region  and  w = 0.25A  (see  Figs. 

4 and  5) . 
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Table  2.  Reactive  load  value  (B.  .)  and  short  circuit  distances 

Load 

(d  (0))  for  a nine  element  aperture  array  using  six  strips 
per  aperture  region,  waveguide  WR-75  (2.159  cm  x 1.206  cm  - 
wall  thickness  0.127  cm),  W = 0.25X,  = 1.55,  and  frequency  • 

7.88  GHz  (see  Fig.  7). 


— 

Aperture 

No. 

0 

= 45® 

0 = 

60® 

®Load 

d^(6)* 

'^Load 

d^(6) 

1 

+0.070 

ml/ 

26.041  mm 

-1.992  mU 

6.392  mm 

2 

-1.56 

7.381 

+0.219 

26.658 

3 

+0.137 

26.319 

-0.457 

11.007 

4 

-0.410 

11.193 

-0.584 

10.513 

6 

-0.658 

10.233 

-0.637 

10.312 

7 

-0.429 

11.118 

-0.578 

10.536 

8 

-0.453 

11.023 

-0.279 

11.722 

9 

-0.749 

9.898 

-0.306 

11.612 

Aperture 

No. 

0 = 75® 

0 = 90 

0 

Load 

d^(0) 

®Load 

d^(0) 

1 

+0.146  ml/ 

26.356  mm 

+0.181  mV 

26.501  mm 

2 

-0.456 

11.011 

-0.468 

10.964 

3 

-0.457 

11.007 

-0.422 

11.146 

4 

-0.575 

10.548 

-0.617 

10.388 

6 

-0.521 

10.756 

-0.617 

10.388 

7 

-0.459 

10.648 

-0.422 

11.146 

8 

-0.207 

12.018 

-0.468 

10.964 

9 

-0.209 

12.009 

+0.181 

26.501 

* Distances  are  measured  from  the  aperture  Interface  or  a multiple 
of  one-half  wavelength  from  the  aperture  plus  the  given  short 
circuit  distance. 


VI.  DISCUSSION  AND  CONCLUSIONS 


It  was  found  that,  for  lossless  aperture  arrays  consisting  of 
eleven  elements  or  less,  the  equivalent  magnetic  current  sheet  should 
be  subdivided  into  three  magnetic  current  strips  for  an  aperture  width 
of  O.lOX,  and  six  magnetic  current  strips  for  an  aperture  width  of 
0.30X.  The  criteria  used  was  the  change  in  magnitude  of  aperture  half- 
space admittance  terms  as  the  number  of  strips  was  changed,  and  whether 

Vis 

or  not  a positive  definite  Real  (Y  ) matrix  resulted.  When  the  number 

of  apertures  increased  beyond  eleven,  the  number  of  required  subdivisions 

was  not  easy  to  predict  for  the  lossless  case.  For  Instance,  for  a fifteen 

element  lossless  aperture  array,  30  subdivisions  were  required  before  a 

Vis 

positive  definite  matrix  Real(Y  ) was  realised.  In  addition,  using  this 
Vis 

[Y  ] matrix  in  an  optimization  subroutine  proved  to  be  costly  in  computer 
(CPU)  time,  due  to  slow  convergence. 

Adding  loss  for  aperture  arrays  consisting  of  eleven  elements  or 
more  appears  to  be  the  best  compromise  between  the  maximum  gain  values 
realized  and  minimum  computer  (CPU)  time  required  for  optimization  (Fig. 

8 is  an  example) . The  use  of  the  number  of  subdivisions  required  for  a 
given  aperture  width  for  an  eleven  element  or  less  lossless  aperture  array 
proved  also  to  be  satisfactory  for  a lossy  aperture  array  with  elements 
ntmberlng  greater  than  eleven. 

From  Figs.  4 and  5 it  was  observed  that  the  beamwldth  was  large 
at  the  low  beam  steering  angles.  This  observation  was  true  both  for  the 
case  of  complex  excitation  (all  elements  fed)  and  for  the  case  of  single 
waveguide  excitation  and  reactive  loading.  One  method  of  producing  a 
narrower  beam  at  low  beam  steering  angles  is  to  place  a dielectric  slab  over 
the  apertures.  The  function  of  the  dielectric  slab  is  to  transform  the 
emanating  electromagnetic  wave  into  a surface  wave  which  results  in  increased 
directivity. 

Villeneuve,  Behnke,  and  Rummer  of  Hughes  Aircraft  [3]  attempted  to 
Increase  the  coverage  of  an  eight  element  aperture  array  with  all  elements 
excited  by  narrowing  the  beamwldth  at  low  beam  steering  angles.  They  used 
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a surface  wave  structure  which  consisted  of  5 layers  of  8-mil  thick 
glass  tape.  A narrower  beam  and  higher  pattern  level  was  realized  in 
the  0 - 0 direction,  while  almost  no  change  resulted  in  the  pattern 
In  the  6 = 90°  direction.  Their  results  were  obtained  experimentally. 

Further  work  should  continue  into  establishing  a theoretical 
explanation  of  the  surface  wave  effects  which  improve  the  directivity 
for  a reactive  loaded  aperture  array  at  low  beam  steering  angles,  and 
to  Investigate  two  dimensional  arrays  similar  to  the  linear  array 

studied  here. 
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PART  TWO 
COMPUTER  PROGRAM 


I . INTRODUCTION 

The  required  subroutines  used  to  generate  the  figures  in  this 
report  are  Included  in  two  previous  reports,  with  the  exception  of 
the  subroutine  which  generates  the  half-space  admittance  matrix  for 
a wide  aperture  aiitenna  array  which  is  described  and  listed  in  this 
report.  Report  [1]  contains  a description  and  listing  of  the  sub- 
routines SICI,  PATHS? , MAXGCV,  MAXGRV,  LINER,  BLOADC,  BLOADR,  FUNCTA, 
FUNCTB,  LINEC,  EIGEN,  and  the  calling  MAIN  program,  while  report  [4] 
contains  a description  and  listing  of  a univariate  optimization  program. 
(Note  that  other  optimization  programs  could  be  used  to  find  the  reactive 
loads  required  for  maximum  gain.) 


II.  HALF-SPACE  ADMITTANCE  MATRIX 

The  subroutine  YHSPA  (N,  Nl,  W,L,T,YHS)  computes  and  stores 
columnwise  in  YHS  the  elements  of  the  admittance  matrix 


4w 


2t  2 

n Lp 


I ^ 


^(iq) (jr) 


(1) 


where 


^(lq)(jr)  4^7t  ~ 


[2S1  Si(u2)  - SKV2)] 


U2  - k 
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Ci  (x) 


Cl  (x) 


cos  V 


dv 


sin  V 


dv 


L Is  the  length  of  the  aperture  In  wavelengths 

separation  distance  between  the 

1th  aperture  q-current  strip  center  and 
the  jth  aperture  r-current  strip  center. 
When  1 = j and  q = r, 


The  following  Is  a list  of  the  abbreviations  used  in  the  subroutine 
calling  statement: 

N is  the  number  of  apertures  In  the  antenna  array. 

N1  is  the  number  of  magnetic  current  strips  per  aperture. 

W Is  the  width  of  the  aperture  (inner  waveguide 
dimension)  in  wavelengths. 

L is  the  length  of  the  aperture  (inner  waveguide 
dimension)  in  wavelengths. 

T is  the  wall  thickness  of  the  waveguide  used  in 
wave lengths . 

YHS  is  the  storage  location  for  the  resulting  half-space 
admittance  matrix. 


DO  loop  10  generates  the  between  the  centers 

of  all  the  magnetic  current  strips  of  the  aperture  array.  DO  loop  11 

evaluates  on  the  distinct  x/.  v terms  needed  to  generate  the 

(iq)(jr)  ij 

terms  - there  are  2*N*N1-N-N1+1  distinct  terms.  DO  loop  15  evaluates 

tl8 

Eq.  (1)  using  the  appropriate  weighting  coefficients  for  each  ^(iq)(jr) 

term.  DO  loops  24  and  25  use  the  upper  triangular  terms  of  the  sym- 

hs 

metric  admittance  matrix  [Y  ] to  generate  the  lower  triangular  terms. 
Minimum  allocations  are  given  by 

COMPLEX  YHS(N*N),  YHS1(2*N*N1-N-N1+1) 

DIMENSION  X(N*N1) 
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PROGRAM  LISTING 


SUBknuTINE  YMSPA(N,Nl,W,LtT, YHSI 
COMPLEX  YHS(<VS») , YHSl  ( 72)  ,U 
DIMENSION  X(42) 

REAL  L,L I 

100  Fl)RNIAT(/3X,'X*//I  JX,6F7.3)  ) 

101  FURMATI/2X,'  Y HALE  SRAC  E * //(  3X  , 5E  l^».7  ) ) 

PI=3.1415Si 

PI2  = 2.*P  I 
ETA=376. 730 
L=L/2. 

Ll=L*2.#PI 

0=2. 

0= (G. .1.  ) 

L2  = Li'»Ll 
N2  = N*M 

M3  = 2*N*NI-N-N  14-1 

N4=M*l 

N5=N2*l 

N6  = 2*M-l 

N7=N-l 

W1  =W/M 

TL=ti*Ll 

TL2=TL’I-TL 

CS=CCS(L I ) 

SN=SIN(L 1 ) 

C l = W/(  PI  ’1‘ETA*L*SN*SN<'N1>!‘N1  ) 

C2=SN*CS 

C3=.5*C0SITL) 

X< 11=0. 

IE IN.EQ. I ) GO  TO  9 

XA=Wl*PI2 

DO  10  1 = 2, N2 

XI  I )=(  I-l  )*XA^(  ( ( I-U/Nl  )*2.*T  )*PI2 
10  CONTINUE 
9 CONTINUE 

WRITE(3,IC0)  (Xm,I=l,N2) 

11  = 1 

DO  11  1=1, N3 
IF(I.hO.l)  GO  TO  12 
IF ( I.GT.N2)  GO  TO  13 
D2=>n  )*X(  I ) 

D=SGRT(D2  ) 

GO  TO  14 

12  D=(U1*2.*PI )/4. 

D2=D*0 

GO  TO  14 

13  I2=(  I-Nb)/N7+2 
XJ=X(  I 1)-X(  12) 

D2=XJ*XJ 

0=SCKT(02) 

I1=I1+N1 

14  SI  = SQRT(  D2«-L2  ) 

S2=SQRT(02-^TL2) 

V1=S1+L1 
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Ul*Di/Vl 

U2*S2+TL 

V2*D2/U2 

CALL  SICI (SD,CD,U) 

CALL  SlCl CSUl fCUl fUl I 
CALL  SICI ISVltCVl  ,V1) 

CALL  SICI  (SU2,CU2,U2» 

CALL  SICI (SV2,CV2,V2» 

YHS2=Cl<'IC2*t  SU2-SV2— J*(SVl-SU  ln-C3*(lJ>*(CUl-C0♦CVl)-CJ2-CV2»-CUl> 
lQ♦CC-CVl  I 

YHS3=C1*(C2*C  CV2-CU2*U*ICV l-CU 1 )I -C^* ( U* ( S U l-SU+SV I ) - SJ2-SV2 )-SU  1+ 
IQ*SC-SVII 

YHSIII »=YHS2-U*YMSi 
IF(I.Eg.N2l  Il=N4 
IFdl.eu.NS)  I1  = N<» 
ll  CONTINUE 

WRITE(3,IC1)  (YHSl(II,l=l.N3» 

DO  15  1=1, N 
1 1=(  I-U+N+l 
YHS(ll»=(C.,0.) 

IF(I.GT.l)  GO  TO  16 
Jl=M 

DO  17  J=l,Nl 
IFU.EQ.NU  GO  TO  18 
YHSm=YHS(  U +J*2  .*YHS1(J1  I 
GO  TO  19 

la  YHS(1I=YHS(1I+N1>*YHS1(  Jl> 

19  Jl=Jl-l 
17  CONTINUE 

GO  TO  15 
16  J1=I*M 
J2=N3-N+I 
DO  20  J=1  ,No 
IFU.GT.NII  GO  TO  21 
YHSim  = YHS(  I ll+J*YHSlIjn 
GO  TO  22 

21  YHSI 1 1»=YHSI 1 1 )+l J-Nl) ♦yHSl(J2» 

J2=J2-N7 

22  J1=J1-1 

20  CONTINUE 
15  CONTINUE 

IFIN.EO.l)  GO  TO  23 
DO  24  I = 2,N 
I 1=1 

DO  25  J=l,N 
J1=N*IAB  sn-j)>i 
YHStm  = YHSUl) 

I 1 = 1 1 + N 
25  CONTINUE 
24  CONTINUE 

23  CONTINUE 
RETURN 
END 
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